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ABSTRACT 


An  alternative  to  the  conventional  rigid  wing  was 
explored  to  see  if  any  advantage  sight  be  obtained  by 
utilizing  slotted  airfoil  sails  for  wings  and  a 
complete  shift  in  emphasis  towards  a lighter  vehicle 
at  the  expense  of  drag  penalties.  Optimum 
combinations  of  size  and  structural  parameters  were 
sought  through  extensive  computer  analysis.  Various 
tradeoff  considerations  for  further  improving  aircraft 
performance  were  considered.  Predictions  of  optimum 
structural  parameters  were  obtained  with  aircraft 
power  requirements  and  the  optimum  flight  velocity. 


4 


TABLE  OF  CONTENTS 


I.  NATURE  OF  THE  PROBLEM 7 

A.  HUMAN  ELEMENT 8 

1.  Power  Versus  Endurance 8 

2.  Krener  Competition 9 

3.  Optimum  Power  and  Velocity 9 

4.  Glider  Philosophy  and  Man-Powered  Aircraft. 

11 

B.  INITIAL  AIRCRAPT  DESIGN 13 

1 . Sail-Ring 14 

2.  Fuselage 14 

3.  Liaiting  Factors 17 

II.  THEORETICAL  DEVELOPMENT 18 

A.  FUSELAGE 18 

1.  Height  Assumptions 18 

2.  Man  and  Propeller  Location 18 

3.  Eapennage  and  Tail-Length  Assuaption s. . . . 19 

4.  Rationale  for  Tube-Spars 20 

B.  HING 22 

1.  Force  Location  and  Identification 22 

2.  Hing-Spar  Coluan  Analysis 22 

3.  Sail  Tensions  and  Allowable  Trailin g-Edge 

Deflections 24 

4.  Location  of  the  Center  of  Gravity 28 

5.  Initial  Height  and  Balance  Caluclations . . 30 


a.  Suaaation  of  Height=Gross  Height  (GH) . 
30 

b.  Hing-Lift (HL)  Equals  the  Gross  Height 


Plus  the  Absolute  value  of  Tail  Lift  (TL) 30 

c.  Balance  of  Aerodynamic  Forces 30 

d.  Height  and  Balance 31 


6.  Spar  Heights  in  Terns  of  Euler's  Buckling 

Equation 31 

7.  Conpressive  Load  Analysis 32 

8.  Spar  Heights  in  Terns  of  Buckling  Criteria 

and  Force  Analysis 33 

9.  Total  Height  and  Wing-Lif t/Unit  Span 

Developnent 34 

10.  Flexibility  of  Wing  Sail  Design 35 

C.  POWER  REQUIRED 35 

1.  Parasite  Drag  Considerations 35 

2.  Induced  Drag  Considerations 38 

D.  SUMMATION 38 

III.  SEARCH  FOB  OPTIMUMS 39 

A.  COMPUTATIONAL  ANALYSIS 39 

B.  COMPUTER  ANALYSIS 39 

1.  "Plot  P"  Optinization 39 

2.  "Tab"  Outputs 40 

3.  Variation  of  Paraneters 40 

C.  TRADE-OFF  ANALYSIS 43 

1.  Two-Man  Crew....... 43 

2.  Magnesiun  versus  Aluniniun 44 

3.  Parasite  Drag  Sensitivity 44 

4.  Sail  Weight  Sensitivity 45 

5.  Taper-Ratio  versus  Pointed  Wingtip 46 

a.  "Weight-Balance"  Conputer  Progran....  47 

D.  AREAS  FOR  FURTHER  CONSIDERATION 48 

1.  Propeller  Optinization 49 

2.  Stability  and  Control 50 

IV.  CONCLUSIONS 51 


So. 





LIST  OF  FIGURES 


1.  Power  Available  vs.  Log  Time 9 

2.  Glider  Vector  Analysis 12 

3.  Ring  Span  and  Sail  Arrangement 15 

4.  Fuselage  Configuration,  Balance  of  Lifting  Forces...  16 

5.  Body  Height  Breakdown 19 

6.  Hing  Spar  Compressive  Load  Analysis 21 

7.  Allowable  Deflections 25 

8.  Selection  of  the  Center  of  Gravity  Location 29 

9.  Parasite  Drag  Coefficient  Cp 37 

10.  Ground  Effect  Considerations 37 

11.  Induced  Drag  Reduction  Polynomial 38 

12.  Graphical  Effect  of  Variation  of  Parameters 41 

13-  Tabulated  Effect  of  Variation  of  Parameters 43 

14.  Taper  Ratio  Effects  on  Drag  Increase 47 

15.  Variation  of  Propeller  Effectiveness  with  Diameter..  50 


6 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


A constant  of  Integration,  collection  of  terms 
c*  angle  of  inclination 
AR  aspect  ratio 
b wing  span 

B constant  of  integration 

Bi  collection  of  terms 

B collection  of  terms 
o 

c wing  chord,  column  end  conditions 

C drag  coefficient 
d 

C,.  induced  drag  coefficient 
ai 

C,  parasite  drag  coefficient 
do 

Cj  jib  sail  chord 

C^lift  coefficient 

wing  lift  coefficient 

Clt  tail  lift  coefficient 

c main  sail  chord 
m 

c root  chord 
o 

Q 

p parasite  drag  coefficient 
d midspan  sail  deflection 
e Oswald's  efficiency  factor 
E modulus  of  elasticity 
f parasite  drag  area 
GW  aircraft  gross  weight 
h spreader  distance 
h/b  height  above  ground/wingspan 
I moment  of  inertia 

K reduction  of  induced  drag  due  to  ground  effect 

K,  deflection  constant  of  proportionality 
a 

tail  area  constant  of  proportionality 

L distance  from  center  of  gravity  to  lifting  force  of  wing 
w 

P power,  pressure 


P power  including  ground  effects 
8 

P inner  compressive  load 

P outer  compressive  load 

o 

P critical  compressive  load  for  column  buckling 
cr 

P interpolating  polynomial 
n 

q dynamic  pressure 
Q inverse  of  slenderness  ratio 
R radius  of  curvature 

outer  restraining  stay 
R2  inner  restraining  stay 

r outer  spar  radius 

0 

r inner  spar  radius 

S sail  area 

SF  safety  factor 

S tail  area 
t 

S wing  area 
w 

(5e  circumferential  stress 

^ sail  weight  per  unit  area 

T tension  from  force  of  sail 

t spar  thickness 

TL  tail  lifting  force 

Tj  jib  sail  tension 

T main  sail  tension 
m 

v velocity 

W weight  of  aircraft 

W body  weight 
b 

inner  spar  weight 
W outer  spar  weight 
WL  wing  lifting  force 

WL'  wing  lifting  force  minus  sail  weight 
X ratio  of  outer  to  inner  wing  spar  radii 
y distance  along  spar 
z deflection  above  span 
Z slenderness  ratio 


I.  KATIJRE  OF  Xfil  PROBLEM 


For  centuries  nan  has  dreaned  of  the  freedom  of  flight 
under  his  own  power.  The  prospect  has  received  much 
attention  and  effort  throughout  the  history  of  mankind.  Yet 
since  the  time  that  Icarus  spread  his  wings  above  the  shores 
of  Crete,  such  endeavors  have  had  little  more  success  than 
that  of  the  Athenian.  This  century,  however,  has  seen  the 
development  and  refinement  of  lighter , stiff er  materials 
coupled  with  an  increased  and  quantified  understanding  of 
aerodynamics.  This  combination  of  increased  awareness  and 
ability  to  solve  the  problem  resulted  in  a string  of 
marginally  successful  man-powered  aircraft.  One  common 
characteristic  of  these  machines  was  a rigid  wing  of  one 
form  or  another.  The  possibility  exists  that  other  wing 
designs  more  suitable  to  the  low-speed-flight  regime  and 
stringent  weight  requirements  of  man-powered  flight  might  be 
feasible. 

While  man  has  been  dreaming  of  flight  for  centuries,  he 
has  been  simultaneously  plying  the  oceans  of  the  world  under 
the  power  of  sail-driven  vessels.  Sails  provided  the  means 
to  drive  ships,  generally  in  winds  of  moderate  force  and 
speed.  Sails  have  also  been  used  successfully  for  centuries 
on  high-powered  windmills.  Modern  sails  are  both  strong  and 
quite  light.  Could  the  technology  of  modern  sails  and 
modern  lightweight  materials  provide  a lighter  flexible  wing 
than  current  rigid  designs?  A Jibsail  and  mainsail 
combination  rotated  to  the  horizontal  provides  a 
configuration  resembling  a slotted  airfoil.  Could  this 
combination  produce  sufficient  lift  and  tolerable  drag  at 
acceptable  power  levels  to  produce  a viable  alternative  to  a 
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rigid  wing?  Is  there  a difference  in  rigid  wing  performance 
and  sail  performance  in  this  low-speed  regime?  These 
intriguing  questions  were  the  subject  of  extensive 
consideration. 


A.  HOMAN  ELEMENT 


1 . Power  Versus  Endurance 


Foremost  in  the  consideration  of  man-powered 
aircraft  is  the  human  being's  ability  to  supply  only  a 
finite  amount  of  power.  This  power  must  be  used  to  support 
the  weight  of  the  pilot  and  his  craft  as  well  as  to  propel 
it  through  the  air.  References  6 and  8 showed  man's  power 
output  ability  as  a function  of  exertion  time  (Figure  1) . It 
clearly  dictates  the  power  for  which  a craft  must  be 
designed  to  operate  for  any  specified  flight  time. 


t .•  ' * M 


Figure  1 -POWER  VS»  LOG  TIME 
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2. 


In  order  to  establish  specific  goals  for  the 
perfornance  of  a nan-  powered  aircraft,  the  requirements  set 
forth  by  Kremer,  a British  industrialist.  They  are  quoted 
in  part  from  Ref.  8. 


A prize  of  L10,0Q0  will  be  awarded  to  an  entrant  from 
any  part  of  the  world  who  first  fulfills  the  conditions. 
The  course  shall  be  a figure  of  eight  embracing  two 
turning  points  which  shall  not  be  less  than  1/2  mile 
apart.  It  must  be  ensured  that  the  machine  is  in 
continuous  flight  over  the  entire  course  and  must  be 
flown  clear  of  and  outside  each  turning  point.  Ground 
clearance  will  be  minimum  of  10  feet  at  start  and  also 
at  the  finish,  both  to  which  are  the  same  point  halfway 
between  the  turning  points.  Between  start  and  finish 
the  ground  clearance  is  unrestricted. 

Based  on  the  Kremer  requirements  and  the  limitations 
of  Figure  1 a designed  flight  requirement  of  .3-. 4 
horsepower  was  selected.  The  craft  had  to  be  maneuverable 

o 

enough  to  negotiate  turns  in  excess  of  180  and  operate,  at 


least  briefly,  at  or  above  ten  feet  ground  clearance. 
Aircraft  design  proceeded  with  these  objectives  and 
limitations  in  mind. 


3.  Optimum  Power  and  Velocity 

The  first  derivative  of  an  equation  set  equal  to 
zero  enables  one  to  obtain  an  expression  for  certain 
optimums  within  that  equation.  In  this  instance  the 
derivative  of  the  power  equation  was  taken  with  respect  to 
velocity  and  set  to  zero.  From  this  an  expression  for  the 
optimum  velocity  was  obtained  and  when  this  was  reintroduced 
into  the  power  equation  an  optimum  power  with  respect  to 
velocity  was  obtained.  This  was  done  as  a first  step  in 
searching  for  the  most  desireable  flight  regime  and  best 
obtainable  power  requirements.  The  development  below 
develops  expressions  for  the  optimum  velocity  and  the  power 
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required  at  the  optima  velocity.  In  the  development  the 
lift  developed  by  an  aircraft  was  assumed  to  just  equal  its 
weight . 


C C ♦ Cl 

d=  do  

T^ARe 


P=power  to  propel  aircraft 
V/*  ir^L  S 


2 3 2 , ^ i/?t  , 

>=<^v  C^Sv/2  =.5^Sv  [C  *C1  / ARe]  *. 


ainiaua  power  occurs  at  dP/dv=0  or: 


3v  C =H  v /[  S TTAR  e (J®v  /2)  ] 


therefore:  P=Pv  S4C  /2 

> do 


. 5ov  =H/sJ3C  ir  ARe  P=4SvC  W/SV3C  “ItARe 

V " d°  do  do 

The  first  derivative  was  solved  and  an  expression  for 
velocity  was  substituted  back  into  the  original  power 
equation  which  was  subsequently  regrouped. 


2H/^S/  (3Cd"JTARe) 


P 4C  H II  2W  1 

= —i? — V—  :» 

3C  TTARe  «DS  (3C  ARe) 
do  ' do 


Sr.- 


(3Fe) 


.25  .25 

V 

AR . 7 * 


AR=b  /S 


4 (2/p)  (f)  (H/b) 


(3TTe) 


V _ (2/p)  (W/b)  '3 

^ .25  .25 

(3He)  (f) 


optiaua  power  and  velocity  teras 

Terms  of  significance  in  both  equations  included  the 

parasite  area  "f"  (where  f*s  times  Cd  ) and  wing-lift  per 

o 

unit  span  (WL/b) . In  both  equations  the  parasite  area 

appeared  to  be  fairly  insensitive  to  change  and  therefore 
not  the  subject  of  concern.  Its  effect  was  examined  in 
computer  analysis  and  receives  attention  subsequently.  The 
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ter*  WL/b  was  seen  to  more  significantly  affect  the 
power-required  equation.  Therefore  the  search  for  the 
configuration  with  the  least  power  became  a search  for  the 
airplane  with  the  lowest  HL/b.  The  most  desirable  locale 
for  the  demonstration  is  at  sea-level  and  therefore  the  air 
density  was  treated  as  a constant.  Oswald's  wing  efficiency 
factor  (e)  was  the  final  variable  to  be  examined  for 
obtaining  optimum  power. 

4 . Glider  Philosophy  and  Man-Powered  Aircraft 

All  of  the  man-powered  aircraft  to  date  have  put 
obvious  emphasis  on  low  drag  and  streamlining  even  at  an 
obvious  cost  in  weight.  Any  knowledgeable  reader  will 
wonder  who  is  wrong,  particularly  when  the  best  gliders  are 
the  fast  streamlined  ones.  To  understand  the  comparison  it 
is  necessary  to  first  better  understand  the  glider  problem: 
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Pigure  2.  GLIDES  VECTOR  ANALYSIS 
Because  of  the  balance  of  forces  and  the  above 
geometry  it  can  be  seen  that  the  glide  slope  or  angle  is 
equal  to  Drag/Lift.  A glider  pilot  simply  wants  to  stretch 
his  glide  as  far  as  he  can  and  consequently  wants  the 
flattest  angle  possible. 


Tying  this  to  the  algebraic  in  simplest  concept: 

D^.5^v  /Wl»Js 

L .5^v  C^S 

In  this  equation  Cd  (parasite  drag  coefficient)  nay 
be  considered  to  be  constant  but  Cl  and  V both  vary  when  the 
pilot  tries  to  fly  at  different  air  speeds.  They  are 


2 

related  since  L=W=.5^Jv  C^S  and  this  is  used  to  eliminate 
and  retain  only  7 in  the  equation. 


2 

D . 5 pv5( Cd 
L~  V 


2 2 
(W/.  5 V S)  \ 
-rrf&<z  ) 


. 5(>v 


. 5pv  SCd  ♦ 


v S 


This  can  be  aininized  by  setting  the  first  derivative  of 
(D/L)  equal  to  zero 


9 vSCd  - 2W 
o 


TTah.5  v S 


2 2 
. 5^v  SCd  =tf/TTAR.5^v  S 

Or  in  words,  the  best  glide  angle  is  at  the  velocity  which 
■akes  the  induced  drag  equal  to  the  parasite  drag.  If  the 
glider  pilot  flies  faster  than  this  the  parasite  drag  goes 
up  wore  than  the  induced  drag  is  reduced.  Slower  flight 
results  in  the  reverse. 

Interestingly  enough,  if  the  glider  pilot  adds 
weight  to  his  airplane  he  can  still  fly  at  the  sane  glide 
angle: 

2 

D/L*Cd/Cl*  (Cd  +C1  /WkR|/Cl  = Cd  /Cl+Cl/TVifitf 
o o 
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If  the  glider  is  flown  at  the  same  lift  coefficient  the  same 
glide  angle  will  result.  With  the  extra  weight  the  craft 
will  fly  a little  faster  to  regain  the  sane  lift 
coefficient. 

If  the  glider  pilot  were  to  improve  his  craft  by 
increasing  the  span  then  that  would  reduce  the  induced  drag 
term.  Extra  advantage  of  this  can  be  taken  by  flying  a 
little  slower  than  previously  to  also  decrease  parasite 
drag. 

Fairings  nay  instead  be  added  to  inprove  the 
airplane  in  order  to  reduce  the  parasite  drag.  The  weight 
of  the  fairings  does  not  penalize  the  glider  but  the  reduced 
drag  is  of  benefit.  Further  benefit  can  be  obtained  by 
flying  a bit  faster  to  reduce  the  induced  drag.  In  sunnary 
drag  and  aspect  ratio  are  important  but  weight  is  not  very 
important. 

Where  is  the  man-powered  airplane  different?  The 
energy  that  a nan  puts  into  sustaining  flight  is  equal  to 
the  airplane  drag  times  the  velocity.  Relating  this  to  the 
glider:  power=(D/L)  tines  V times  W.  The  D/L  indicates 
some  of  the  characteristics  of  the  glider  but  in  addition 
must  consider  the  W and  V terns.  In  the  glider,  additional 
weight  resulted  in  slightly  faster  flight  to  regain  the  same 
D/L  , but  in  this  instance  both  W and  V will  increase  for  a 
similar  situation.  On  the  other  hand  weight  reductions  mean 
both  W and  V would  decrease  even  with  reduced  velocity 
sufficient  to  maintain  the  s ame  glide  angle.  Apparently 
this  kind  or  argument  has  not  been  recognized  by  some  of  the 
other  competitors  for  the  Kremer  prize  because  they  all 
seem  to  be  trying  to  design  gliders  with  propellers.  Weight 
is  sacrificed  in  order  to  reduce  parasite  drag.  The 
opposite  is  true  of  this  analysis. 


B.  INITIAL  AIRCRAFT  DESIGN 
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1 . Sail-Wing 

For  reasons  detailed  in  a subsequent  section  the 
following  wing-spar  and  sail  arrangements  were  selected 
(Figure  2) . The  semi-span  of  each  wing-spar  consisted  of 
two  hollow  cylinders  of  unequal  radius  but  of  equal  length. 
The  inboard  cylinder  was  fixed  at  the  fuselage,  stabilized 
at  its  extremity  by  a wire  stay,  and  joined  to  the  outboard 
cylinder  in  a pin-ended  fashion.  The  outboard  cylinder  was 
restrained  at  the  wingtip  by  a second  wire  stay  which  joined 
the  fuselage  coincidentally  with  the  inner  fuselage 
restraining  wire.  Two  vertical  spreaders  with  stays  were 
positioned  at  mid-span  so  as  to  provide  structural  integrity 
to  the  wing-spar  segments.  A third  vertical  spreader  was 
positioned  on  top  of  the  fuselage  to  support  the  wing 
structure  under  static  conditions.  The  sail  arrangement  was 
chosen  to  resemble  a genoa  jib  which  overlaps  50%  of  the 
mainsail  and  a mainsail  combination.  Hhen  this  combination 
is  rotated  from  the  vertical  (sailing  orientation)  to  the 
horizontal,  one  wing  of  a slotted  airfoil  is  obtained.  The 
initial  analysis  considered  a wing  with  a pointed  wingtip 
and  subsequent  analysis  examined  the  trade-offs  with  a 
finite  chord  tip. 

2 . Fuselage 

The  fuselage  design  consisted  of  a single  hollow 
tube  referred  to  as  the  fuselage  spar  (Figure  4) . The 
propeller  was  placed  at  the  forward  extremity  of  the 
fuselage  with  the  man  arbitrarily  positioned  two  feet 
further  aft.  Shortly  behind  the  propeller  the  leading  edge 
of  the  jibsail 
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attached  to  the  fuselage  spar.  A specified  distance  aft  of 
the  nan  the  wing-spar  attached  to  the  fuselage-spar. 
Directly  above  the  wing-spar  junction  with  the  fuselage  the 
static  support  spreader  was  located.  Enpennage  positioning 
was  deternined  by  aircraft  weight  and  balance  criteria.  To 
be  as  light  as  possible  the  size  and  thickness  of  these 
aenbers  were  obliged  to  be  optimized  consistent  with 
structural  and  safety  limitations.  These  limitations, 
however,  were  not  the  only  constraining  elements. 

3 • Limiting  Factors 

Primary  structural  criteria  for  spar  design  were 
bending  and  Euler  buckling.  Therefore  for  a given  strength, 
the  spar  grew  lighter  as  the  radiis  was  increased  and 
thickness  reduced.  However  to  prevent  local  crippling  and 
practical  handling,  a minimum  wall  thickness  of  .030"  was 
chosen.  Lift  coefficient  upper  limits  were  not  definitely 
specified,  but  attempts  were  made  to  avoid  exceeding  values 
of  1.5.  Modern  hang-glider  wings  have  achieved  Cl's  greater 
than  2.0  in  tests  and  Cl-stall  in  the  region  of  1.6 
according  to  Ref.  9.  Therefore  Cl=1.5  was  determined  to 
leave  sufficient  margin  for  control,  gusts,  etc.  Last  and 
most  important,  was  the  power-reguired  for  the  designed 
craft.  The  intention  was  to  produce  a prototype  at  least 
capable  of  traversing  Kramer's  course  and  preferably  able  to 
sustain  flight  for  approximately  twenty  minutes  in  duration. 
Consequently  efforts  were  directed  at  achieving 
power-requirements  in  the  .3-. 4 horsepower  region  as 
suggested  by  Figure  1.  These  criteria  served  as  the 
principal  limiting  values  in  the  minimizations  and  trade-off 
considerations  that  followed. 
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II . theoretical  development 


A.  FUSELAGE 


1 . Height  Assun ptions 

In  the  hope  that  this  design  night  be  utilized  by 
■ore  than  the  narrow  spectrum  of  light-weight 
"champion-athletes",  the  weight  of  the  man  was  assumed  to  be 
180  pounds.  Propeller  and  bicycle  apparatus  were  each 
assuned  to  be  ten  pounds.  Empennage  weights  were 
considered  to  be  approximately  20*  of  their  respective  wing 
counterparts  while  the  fuselage  spar  was  initially  assuned 
to  be  about  24  pounds.  The  "body  weight"  breakdown  for  an 
initial  guess  was  selected  at  240  pounds: 


object 

weight 

■ an 

180 

propeller 

10 

bicycle  apparatus 

10 

enpennage 

16 

fuselage  spar 

24 

TOTAL 

240 

.2(wt  +wt  ) 

spars  sail 


Figure  5.  BODY  HEIGHT  BREAKDOWN 
2.  Han  and  Propeller  Location 

Three  options  were  considered  for  the  positioning  of 
the  nan  and  propeller.  First,  the  conventional  arrangement 
with  the  propeller  at  the  extreme  forward  position  and  the 
nan  somewhat  further  aft.  Second, 
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the  aan  in  the  forward  region  of  the  aircraft  but  the 
propeller  at  the  extreme  aft  position  of  the  fuselage. 
Third,  the  man  in  the  forward  region  of  the  fuselage  with 
the  propeller  immediately  to  his  rear.  The  second 
configuration  was  abandoned  due  to  the  complexities  of  a 
propeller  amidst  the  empennage  structure  and  also  due  to 
anticipated  weight  penalties  related  to  the  extended  drive 
mechanism.  Although  attractive  from  a drag  standpoint  (i.e. 
pilot  forward  of  the  slipstream)  the  third  option  was 
discarded  because  of  the  pronounced  center-of-gravity 
effects  and  weight  penalties  of  a man-first  arrangement. 
Also,  the  structural  difficulties  associated  with  a rotating 
propeller  on  a bending  beam  were  considered  undesir eable. 
The  first  option  was  chosen  for  several  reasons.  At  the  low 
speeds  anticipated  for  the  flight  regime  drag  was  not 
considered  to  be  a significant  factor.  This  was  borne  out 
in  mathematical  analysis  of  the  power  equation  and 
subsequent  experimental  results.  In  addition,  the  man  aft 
of  the  propeller  provided  for  a shorter  fuselage  due  to 
center-of-gravity  considerations.  This  meant  a weight 
reduction  which  in  turn  resulted  in  a power  reduction.  A 
man  close  to  the  center-of-gravity  also  meant  a reduction  in 
the  bending  stresses  which  enabled  the  employment  of  a 
lighter  fuselage  spar  (thinner  and/or  smaller  in  radius) . 
While  unconventional  arrangements  of  man  and  propeller  were 
intriguing,  there  appeared  to  be  valid  reasons  for  what  is 
considered  as  the  •,conventional,,  arrangement. 

3.  Empennage  and  Tail-Length  Assumptions 

Por  ease  of  calculations  and  flexibility  the  size 
and  weight  of  the  empennage  were  considered  to  be  20%  of 
their  respective  wing  counterparts.  This  was  a good  first 
approximation  for  the  optimization  analysis  which  kept  the 
number  of  parameters  manageable.  Tail  length  considerations 
were  based  on  the  illustration  in  Pigure  4. 

A balance  of  force  moments  about  the 
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W1 


Eqn  (1) 


center-of-gravity  yielded: 

Cl  qS  (L  )=C1  qS  (L  ) 

WWW  t t t 


Cl  S =C1  S L 
t t w 


V 

t 


V=  1 Qv  (Cl  S -Cl  S ) *‘  P»  Cl  S (1-L  ) 
— \ w w t t — \ w w — r 

2 2 L 


Eqn  (2) 


For  reasons  described  subsequently  "L  M was  assumed 

w 

equivalent  to  c /6  and  "L  '•  was  assumed  to  be  equivalent  to 
o t 

1.5  times  the  wing-sail  root  chord.  This  was  deemed 

adequate  to  provide  a sufficient  tail  volume  for  aircraft 
stability.  Equation  (2)  then  became: 


W= 1 pv  s Cl  1-c  /6 
— \ w w —o — 


1 . 5c 


U lov  S Cl  (1-  1 ) qS  Cl  (.89) 

1 — \ w w 


Eqn  (3) 


2 ' 6 (1.5) 

Alternately,  this  means  that  the  lift  provided  by  the  wing 
(WL)  must  be  12. 5X  greater  than  the  aircraft  gross  weight. 


4.  Baticnale  for  Tube-Spars 

Hollow  tubes  were  selected  for  both  the  fuselage  and 
the  wing  spars  for  several  reasons.  Three  major  structural 
requirements  were  considered.  They  were:  bending  moments, 
buckling  . forces  and  torsion.  Hollow  tubes  were  not 
considered  inherently  superior  to  other  possible  spar  shapes 
but  in  the  fuselage  spar  it  was  the  prospect  of  torsion 
between  the  propeller,  wings  and  tail  section  that  prompted 
a hollow  tube.  In  the  wingspars  torsion  was  a 
consideration;  however  buckling  loads  were  considered  the 
principal  structural  consideration  here.  Hollow  tubes  were 
desireable  for  the  streamlined  effects  of  the  leading  edge 
in  low-speed  flight  regimes.  Overall,  the  choice  of  hollow 
tubes  provided  a versatile  structural  member  for  anticipated 
load  conditions. 
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B 


WING 


1 . Force  Location  and  Identification 

Wing  spars  were  treated  as  columns  for  purposes  of 
structural  analysis.  A schematic  of  compressive  loading  is 
shown  in  Pigure  6.  Compressive  forces  were  observed  in  the 
outer  wing  spar  due  to  the  tension  force  from  the  jib  and 
mainsail  and  from  the  outer  restraining  stay  (R^).  Inner 

spar  compressive  loads  were  due  to  the  total  outer-spar 

compressive  load  plus  a component  from  the  inner  restraining 

stay  (B  ) . 

2 

2.  Wing-Spar  Column  Analysis 

If  the  wingspar  were  to  be  subjected  to  a vertical 

distributed  load  from  the  mainsail  as  well  as  a compressive 

loading,  then  it  would  reguire  analysis  as  a beam-column. 

This  approach  was  not  pursued  for  several  reasons.  First, 

as  the  subseguent  tension  analysis  demonstrated,  tensions  in 

the  sail  were  regarded  as  acting  primarily  in  the  spanwise 

direction.  It  was  also  intended  that  the  mainsail  was  to 

carry  no  vertical  loads  in  the  coincident  region  of  jib  and 

mainsail.  Finally,  if  some  vertical  loading  were  to  result 

it  could  be  disregarded  in  all  probability  according  Ref  2: 

Since  the  critical  load  for  the  beam  column  is  the  same 
as  the  column  critical  load,  it  is  apparent  that  the 
axial  load  which  causes  instability  is  independent  of 
the  lateral  loads.  This  is  true  of  all  beam  columns  if 
none  of  the  stresses  exceeds  the  proportional  limit. 

All  angles  in  Figure  6 were  considered  sufficiently 

small  for  small-angle  approximations  in  analysis.  The 

weight  of  the  wing  spars  acting  downward  provided  some 

assistance  to  the  restraining  stays  in  balancing  the  forces 

from  the  sail  tension.  Spar  weighty  then  aided  in  reducing 
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compressive  forces  within  the  spars  due  to  the  restraining 
stays.  This  fact  further  assisted  in  weight  reduction. 
From  statics  and  Figure  6: 
at  the  wingtip 


ffFy:  T(<*+  h ) B f h |H 

b/2  ~ [b/q« 


Eqn  (4) 


g-px:  p =t*r 

o 1 


Egn  (5) 


at  midwing 


: M_LL«2LiAl‘V\ 

lb/4  J lb/4|  4 


Eqn  (6) 


£Fx:  P =P  +H 

i o 2 


Eqn  (7) 


for  the  critical  compressive  loading  of  a column  as  follows: 


Eqn  (8)  P =cTTEI 

cr  

2 

L 


where  c=end  conditions 


for  a hollow  tube: 

♦ « 

Eqn  (9)  I=rr(r  -r  ) 

~ o i 
4 

2 2 2 2 
(r  ♦r.  ) (r  -r.  ) 

o 1 o 1 


=1.0  pin  ended. free 
=2.04  pin  ended,  fixed 
E=modulus  of  elasticity 

I=moment  of  inertia 

L=column  length 

r =outer  radius 

0 

r =inner  radius 

1 


2 2 

(r  *r  ) (r  +r  ) (r  -r  ) t=wall  thickness 
o i o i o i 

2 2 

for  tubes  with  large  radii  compared  to  thickness  (r  ♦r  ) 

o i 

2 

is  approximately  equal  to  2r  , where  (r  +r  ) approximately 

o o i 

equals  2r  and  (r  -r  ) =t.  these  approximations  into 

o o i 

equation  (9)  yielded: 

2 3 

1=712  (r  ) 2r  t=  r t tt  Eqn  (10) 

— o o o 
4 

For  an  outer  radius  of  3.0  inches  and  a thickness  of  .03 
inches  this  approximation  compares  as  follows: 

1=  1(3.0  -2.97  ) =2.507  in  Eqn  (9)4 
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Eqn  (10) 

Substituting  equation  (10) 


I=ir(3. 0)  (.03)  =2.54  5 in 

(a  difference  of  about  1.49X). 

into  equation  (8)^  yields 

Eqn  (11)  P =cTT  Er  t where: 

cr  a 

2 

(b/«)  r =outer  radius 

o 

b/4=coluan  length 

This  is  the  maximum  compressive  load  that  a column 
of  given  radius,  thickness,  length,  specific  end-conditions 
and  type  aaterial  can  withstand. 

3.  Sail  Tensions  and  Allowable  Trailing- Edge 

Deflections 


Reference  7 discussed  circuaf erential  stresses 
experienced  in  a thin-walled  cylinder  and  was  directly 
related  to  the  sail  tension  equation  development  illustrated 
in  Pigure  7. 

where 


Eqn  (12) 


6=Pr 
® t 


=wall  tension 
unit  area 

Oot  =wall  tension 

w 

unit  depth 


<$3=circuaf  erential 
stress  in  a thin-walled 

cylinder 
P=pressure 
r=cylxnder  radius 
t =wall  thickness 

w 

R=radius  of  curvature 


Sail  tensions  were  assumed  to  exist  along  the  spanwise 
elements  of  a wing  instead  of  in  the  chordwise  direction. 
Consequently  the  force  of  the  sail  could  be  regarded  as  a 
series  of  tension  strips  running  from  root  to  tip  with 
virtually  nc  fore-and-aft  stress  components.  If  the 
curvature  of  the  trailing  edge  of  the  jib  sail  were  regarded 
as  a segment  of  a very 
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large  circle,  then  the  spanwise  tension/unit  chord  (T/c=t) 
was  equivalent  to  (t) . Tension/unit  chord  was  assumed  to  be 
a constant  along  each  spanwise  strip.  In  the  small  angle 
approximation  the  reciprocal  of  the  sail  radius  of  curvature 
was  known  to  be  equivalent  to  the  second  derivative  of  the 
deflection  (z)  with  respect  to  spanwise  position  (7) . These 
conditions  modified  Equation  (12)  to  the  following: 


T/c=PF 

Por  the  jib  sail: 
2 


1 d z Pc=-Pc  ( 1 — 2 y ) 
. = j - 

H dy  T b 

J 


Eqn  (13 
Eqn  (14) 


where  c^=c(1-2y/b) 


The  negative  sign  originated  from  the  radius  of  curvature 
orientation.  Integrating  equation (14)  twice  and  solving  for 
the  constants  of  integration  produced  the  following  results: 


dz=Pc  (by-y  )♦* 
dy  T 2 

1 2 , 

z=Pc  (by  -y  ) +iy*B 

n 1 1 

with  end  conditions: 


y=0,  z=0 
y=b,  z=0 


B=0 

A=Pc  b 


6T  ' 


Equation  (14)  then  became  an  equation  for  the  spanwise  sail 
deflection: 


9 2 2 

z Pc  (2y  -3by  +b  y ) 

j 


Eqn  (15) 


Solving  Equation  (15)  for  the  spanwise  deflection  at  midspan 
(y=b/4)  resulted  in 
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Eqn  (16) 


z 


3 2 2 

Pc  (2y  -3by  *b  y) 
64? 


j 


While  equation  (16)  provided  a general  expression  for 

■id-span  sail  deflection,  the  aaxiaca  value  (d)  that  the 

sail  would  be  allowed  to  experience  regained  to  be 

determined.  It  was  assumed  that  aid-span  deflection  was 

directly  proportional  to  the  product  of  chord  and  a constant 

of  proportionality  (K  ) . The  equation  indicated  that  the 

d 


tension  in  the  sail  and  the  consequent  design  stay  tensions 


and  spar  compressions  increased  linearly  as  the  sail  was 
drawn  tighter  and  aade  flatter.  This  meant  that  there  was  a 
tradeoff  between  baggy  sails  and  high  loads.  Judgement  was 
exercised  to  allow  the  deflection  to  be  limited  to  5%  of  the 
root  chord. 


d=K  c 
d 

it  aidspan  d=K  c /2  Eqn(18)  . 
d 0 


Eqn  (17) 

For  a aidspan  deflection  of  5* 


K became  .10  Substituting  Equation  (18)  into  Equation  (16) 


resulted  in  an  expression  for  the  required  sail  stay  tension 

2 2 

T Pc  b Pb  c Eqn  (19) 


1* 


* 


64K  c /2  32K  c 
do  do 


But  pressure  (P)  was  equivalent  to  the  total  lifting  force 
(VL)  divided  by  the  total  wing  area(S).  Referring  to  Pigure 
6 it  can  be  seen  that  the  weight  of  the  sails  acts  opposite 
to  the  lifting  force.  Since  the  sail  supported  its  own 
weight  plus  the  weight  of  the  entire  craft  it  was  assuaed 
that  the  weight  of  the  sail  should  not  affect  Equation  (19) 
and  should  therefore  be  deducted  froa  the  effect  of  the 
wing-lift: 

2 2 

T (WL-GS)b  c WL*b  c Eqn  (20) 

•»  j« $■ 

32K  S c 32K  Sc 
d o d o 
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For  all  considerations  the  area  of  the  main  sail  directly 
beneath  the  jib  was  assumed  to  provide  no  lift.  It  was 
assumed  that  this  region  merely  provided  the  slotted 
condition  between  the  main  and  jib  sails.  Consequently  for 
structural  load  considerations  the  main  sail  chord  was 
measured  from  the  trailing  edge  of  the  jib  chord  to  the 
trailing  edge  of  the  main  sail  chord.  Bearing  this  fact  in 
mind  a similar  analysis  of  mainsail  tensions  lead  to: 

T HL'b  c Eqn  (21) 

m* » 

32K  Sc 
d o 

The  jib  tension (T  ) and  main  tension  (T  ) combined  to  equal 

j ■ 

the  total  tension  (T)  acting  along  the  wing  spar  because 

c *c  =c  . Figure  7 shows  a spanwise  tension  balance  acting 
m j o 

at  one  wingtip.  When  these  forces  and  angles  were  summed, 

the  result  was  a compressive  force  acting  aloag  the  wing 
spar.  Fore-and-aft  force  components  were  equal  and 
opposite. 

4.  Location  of  the  Ceater  of  Gravity 

The  mean  aerodynamic  chord  of  a triangular  wing  is 

located  one-third  the  distance  of  the  seni-span  (b/6)  from 

the  fuselage.  Therefore  c = 2c  /3.  The  wing  could  also  be 

ave  o 

regarded  (Figure  8)  to  be  a rectangular  wing  of  the  same 
span  and  two-thirds  the  value  of  the  root  chord.  Stability 
considerations  generally  select  the  center-of-gra vity 
location  close  to  the  quarter  chord  position.  In  this 
instance  the  quarter  chord  position  of  the  equivalent 
rectangular  wing  can  be  seen  to  intersect  the  fuselage  at  a 
value  one-third  of  the  wing  root  chord. 


28 


Jm 


29 


5.  Initial  Weight  and  Balance  Caluclat ions 

Several  criteria  were  followed  for  initial  weight 
and  balance  analysis: 

a.  Summation  of  Weight=Gross  Weight  (GW) 

Gross  weight  included  the  initial  weight 
estimates  of: 

GW=W  = W ♦ (W  ♦ S)  (UK  ) Eqn  (22) 

body  spars  t 

where  K =factor  for  the  tail  structure  and  was  assumed  to  be 
t 

0.200. 

b.  Wing- Lift (WL)  Equals  the  Gross  Weight  Plus  the 
Absolute  value  of  Tail  Lift  (TL) 

The  wing-lift  must  support  not  only  the  aircraft 

gross  weight  but  also  an  additional  amount  equal  to  the  tail 

lift  loss  as  discussed  in  the  Empennage  Assumptions  section. 

=0=  qS  Cl  L -qS  Cl  L 
eg  www  t t t 

S Cl  L =S  Cl  L Eqn  (23) 

w w w t t t 

This  substituted  into  Equation  (22)  yields  an  expression  for 

weight  and  also^wing-lif t of  the  aircraft: 

W lov  S Cl  ( 1 — L ) 

» -\  w w — w 

2 L where 

t 

2 

WL  W 1r>v  S Cl  Eqn  (24) 

= *— \ w w 

(1-L  /L  ) 2 

w t 

c.  Balance  of  Aerodynamic  Forces 

Aerodynamic  forces  must  be  balanced  about  the 
center-of-gravity  as  illustrated  in  Figure  4 Thus 

WL  (L  ) ®TL  (L  ) Eqn  (25) 

w t 
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Due  to  the  location  of  the  center-of-gravity  at  one-third  of 
the  root  chord  aft  of  the  leading  edge  and  the  wing-lift 

acting  at  the  aid-point  of  the  root  chord,  the  value  of  "L  " 

v 

was  found  to  be  one-sixth  of  the  root  chord.  "L  " was 
assuaed  to  be  three-halves  of  the  root  chord  value. 

d.  Height  and  Balance 

After  the  balance  of  aerodynaaic  forces  was 
obtained,  weight  shifts  were  necessary  to  insure  that  the 
weights  actually  balanced  at  the  center-of-gravity . The 
geoaetry  of  the  aircraft  was  fixed  with  the  exception  of 
pilot  location.  Therefore  all  that  reaained  was  to  position 
the  pilot  appropriately  to  bring  the  C.G.  to  the  desired 
location.  Care  was  required  to  insure  that  the  propeller 
was  further  forward  than  the  leading  edge  of  the  jib  sail 
stay. 

6 . Spar  Heights  in  Ter as  of  Euler 's  Buckling  Equation 

A column  subjected  to  compressive  loads  will  fail 

when  a critical  load  is  reached  or  exceeded.  Equation  (11), 

was  utilized  for  the  development  of  wing  spar  weight 

expressions.  Spar  weight  is  equivalent  to  the  material 

density  tines  the  volune  of  the  structure.  Spars  are 

hencefotth  considered  as  "inner"  spars  and  "outer"  spars. 

Por  example,  the  weight  of  the  outer  spars  is  the  weight  of 

both  outer  wing  spar  segments.  Thus 

spar  weight=density  tines  volume*  2 rtb/2  Eqn(26) 

for  an  outer  radius  of  3.0  inches  and  a thickness  of  .03 

inches  this  approximation  compared  as  follows: 

A = TT(r  -r.  )=  (3.0  -2.97  )=.5626 

A =2'*tt*  2 3.  OTT.  03*.  5655 
2 

(a  difference  of  about  0.51X).  Equation  (11)  may  be 
rearranged  to  include  spar  weight  terms  as  follows: 
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3 3 2 2 

P c TTEr  t c^Er  (2^rt^b/2) 

Cr  2 2 

(b/4)  (b/4)  ^b 

2 2 2 2 

P c-n-Er  W 16c  "n'E  (r/b)  W Eqn  (27) 

cr* -spar*  — spar 

(b/4)  ^b  ^b 

while  the  distance  in  the  numerator  is  a semi-span  (b/2)  for 
weight  calculation,  the  denominator  remains  the  length  of 
each  individual  column  (b/4)  for  buckling  considerations. 
The  weight  of  the  wing  spar  is  equivalent  to  the  weights  of 
the  outer  and  inner  wing  spar  sections.  The  (r/b)  term  was 
called  the  "slenderness  ratio"  and  was  one  of  the  terms 
utilized  in  subsequent  minimization  anal/ses. 

7.  Compressive  Load  analysis 


Equations (4-7)  state  compressive  loads  in  terms  of 

aircraft  geometry  (Figure  6) . The  angle  "alpha"  was  assumed 

to  be  equivalent  to  2d/b/4.  From  the  "allowable  deflection" 

discussion  d=K  c.  At  midspan  the  chord  value  was  half  of 
d 

the  root  chord  and  therefore  d=K  c /2.  Looking  at  sail 

d o 

geometry  S=bc  /2  and  therefore  d=K  S/b.  Consequently 
o d 


=8d/b=8K  S/b 
d 


Eqn  (28) 


Solving  Equation  (4)  for  the 


tension  force  in  the  outer  stay  became 


2 


R WL'  b 
1 * 

32K  S 
d 

[2K  S 1 
~d  ♦— 

l bh  2 J 

fob 

16h 

Eqn 

(29) 

from 

Equation  (5) 

the 

compressive 

loading  in  outer 

wing 

spars 

became 

2 

P WL' b ! 

o=  1 

32K  S 

A 

2KS  3 
["bh  2 

1 Wob 
1 16  h 

Eqn 

(30) 

Similarly  from  Equations  (6)  and  (7)  the  following  analysis 
was  obtained  for  the  inner  wing  spars: 
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HL'b  pK^S  1]  b(2Bo+8i) 

2~32K  S [ bh  *2\  ”l6h 
d 

2 

P.  HL'b  pKjS  2]  b(3Ho*Hi) 

1_32fc  S [ bh  * ] 16h 
d 


2qn  (31) 


Eqn  (32) 


8.  Spar  Heights  in  Terms  of  Buckling  Criteria  and  Force 
Analysis 

Euler's  buckling  equation  provided  an  expression  for 
the  critical  buckling  load  that  Bay  be  applied  to  a column. 
It  was  then  aodified  to  reflect  the  weights  of  inner  and 
outer  wing  spars  (Equation  27) . This  set  of  equations  was 
set  equal  to  the  compressive  loads  calculated  from  the 
aircraft  geometry.  Spar  weights  were  calculated  from  these 
equations  after  the  inclusion  of  a safety  factor  (SF)  term, 
a universally  accepted  engineering  practice.  Equation  (27) 
was  thus  modified  to  include  a safety  factor  term  in  the 
denominator.  When  the  safety  factor  exceeded  its  minimum 
value  of  1.0,  the  resultant  decrease  in  tolerable  critical 
loading  provided  a margin  of  error.  Thus,  for  the  outer 


wing  spars: 


: P =P  or 

cr  o 

2 2 2 
1 6c  TTE  (r  /b)  Ho  HL'b  [2K^S  3]  Hob 

SF  b°  32K  S [bh  *2]  iTh” 

d 1 


Solving  for  the  weight  of  the  outer  spars  (Ho)  yielded: 


HL'  1 3bh 
Ho  l 4K  S 


2 2 2 . 
14-16  Z Eh 


Eqn  (33) 


SFpb 


where  Z *r  /b 
o o 


For  simplification  of  calculations  let: 


2 2 2 2 

A bh  B 16  TTZ  Eh  B 16TTEZ  Eh 

o= © i=  i 

2 2 

K^S  SP^b  SF^b 

then  Equation  (33)  became 


R RL* (1 +3 A/4) 

o=— 

(1  + B ) 
o 

Solving  for  R yielded 
i 


H RL'  [-  1 + 2B  + AB  - 5 A/4*] 
i=  o o J 


(1*2B.)  (HB  ) 
i o 


Eqn  (34) 


Eqn  (35) 


9 . Total  Rejqht  and  Rinq-Lif t/Unit  Span  Developme nt 

It  was  determined  from  Equation  (2)  that  aircraft 

weight  was  equivalent  to  RL (1-L  /L  ).  Aircraft  gross  weight 

w t 

was  assumed  to  be  a body  weight  (R  ) , which  was  estimated  in 

b 

a previous  section,  plus  the  weight  of  the  wing  spars  and 
sail,  plus  the  empennage  (assumed  to  weigh  20%  of  their 
respective  wing  components) . Thus 


R = RL(1-L  /L  )=R  M1+K  ) (R  ♦ S)  Eqn  (36) 

w t b t spars 


The  combined  weight  of  the  wing  spars  became 


H =R  + R RL'  [-A/2  + 2(B  + B )+A(B  + 3B  /2)]  Eqn  (37) 

sp  o i= © — i © i 

(1*2b.)  ( 1 + b ) 
l o 

Rhen  this  was  substituted  into  Equation  (36)  an  expression 
for  wing-lift/unit  span  was  obtained. 


s ' 

fR  1+K 

h A 4- 

- 

A/2+2  (B  +B.)+A(B  +3B  /2fl 

1+K 

!)▼  X 

<5s 

(1+2B  ) ( 1 ♦ B ) J 

b 

[(1-L 

/M 

b 

-A/2  + 2 (B  ♦ B )+A(B  +3B  ./2) 

l (HK  )' J 

(1*2B  ) (1  + B*)  1 

Eqn  (38) 


Equation  (38)  represented  the  dominant  term  in  the  optimum 
power  equation  developed  in  a previous  section. 
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Considerable  attention  was  devoted  towards  minimizing  its 
value. 


10.  Flexibility  of  Wing  Sail  Design 

Two  remaining  terms  in  the  optimum  power  equation 
that  could  be  varied  were  the  parasite  drag  area  (f)  and  the 
wing  efficiency  factor  (e) . Both  terms  were  examined  for 
optimization  possibilities.  Aircraft  wing  geometry  might  be 
varied  to  improve  the  wing  efficiency  factor,  referred  to  as 
Oswald's  efficiency  factor.  One  particularly  attractive 
feature  of  Equation  (38)  is  that  it  may  be  optimized 
independently  from  wing  geometry.  Conversely  improvements 
in  wing  efficiency  need  not  alter  wing-lift/unit  span 
calculations  so  long  as  wing  area  remained  fixed. 

C.  POWER  REQUIRED 

Re-examining  the  power  equation  discussed  in  a previous 
section  for  drag  terms  evolved  the  following: 

P*DV1j>v  S(Cd  * Cl  ) Eqn  (39) 

2 P Trise 

where  the  first  term  in  parenthesis  represented  a parasite 
drag  coefficient  and  the  second  term  an  induced  drag. 
Certain  assumptions  as  to  their  nature  preceeded  an 
examination  of  their  minimization. 

1 • Parasite  Drag  considerations 

An  equivalent  parasite  area  (f)  was  obtained  from 

the  product  of  Cd  and  sail  area  S.  This  term  was  further 

P 

broken  down  to  a parasite  drag  for  the  man  and  fuselage 

(referred  to  as  "f  ")  and  a term  relating  to  the  drag  of 

man 

the  wing  spars.  Von  Mises  [Ref.  7]  mentions  a coefficient 


of  parasite  drag  Mc  , which  is  based  on  the  frontal  area  of 

P 

the  body.  Figure  9 [Von  Hises  1959]  presents  a region  of 
expected  parasite  drag  coefficients  for  various  geometries. 
The  anticipated  Reynolds  number  for  man-powered-flight  is  on 

4 

the  order  of  5.5x10  . 


Figure  9-PARASITE  DRAG  COEFFICIENT  Cp 

The  anticipated  Reynolds  number  is  on  the  lower 

extremity  of  this  figure  and  suggests  an  increasing  parasite 

drag  coefficient  value.  However,  a shape  more  streamlined 

than  the  ellipse  was  anticipated  due  to  wrapping  the  sail 

about  the  wing-spar  and  reaffixing  it  to  the  mainsail  some 

distance  aft  of  the  spar.  A value  of  c =.25  was  arbitrarily 

P 

selected  based  on  the  above  factors.  A product  of  wing  spar 
diameter  times  the  wingspan  yielded  the  frontal  area 
mentioned  earlier.  The  parasite  area  thus  became 

f=f  *c  (2r  *2r  } b/2  Egn  (40) 

man  p o i 


Figure  10  - GROUND  EFFECT  CONS I DEBIT 13 NS 
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2.  Induced  Dr ag  Considerations 


Induced  drag  (Cd  ) is  equivalent  to  the  square  of 

the  lift  coefficient  divided  by  the  product  of  pi  and  the 

aspect  ratio.  Figure  10  shows  a reduction  in  the  induced 
drag  and  thrust  required  due  to  a phenomenon  known  as 
"ground-effect.”  Since  the  majority  of  the  man-powered 
flight  regime  was  expected  to  occur  in  close  proximity  with 
the  earth,  ground-effect  provided  an  extremely  useful  drag 
reduction.  With  ground-effect  the  induced  drag  equation 
became 

2 

Cd.  KCl  Eqn  (41) 

"1TAR 

where  K represents  a reduction  due  to  ground-effect.  Osing 
the  curve  of  Figure  10  and  numerical  approximation 

2 3 4 

SDEC  K'  A A A A 

h/b  C (1-X) 
di 


0.1 

.48 

.52 

.19 

o 

• 

|SJ 

.29 

.71 

.10 

-.09 

.05 

0.3 

.19 

.81 

-.04 

-.02 

.06 

.03 

0.4 

.13 

.87 

.05 

-.01 

0.5 

.08 

.92 

INDUCED  DRAG  REDUCTION  FORWARD  INTERPOLATING  POLYNOMIAL 

Figure  11 
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techniques  discussed  in  Gerald  [ Ref . 1 ] a second  order 
polynomial  fit  was  obtained  for  induced  drag  reduction. 
Figure  11  illustrates  the  data  employed  for  the  drag 
reduction  equation  as  a function  of  altitude/wing  span 
P (x)  =K=fo*5£f o*s  (s-1)A  fo  where  s=x- . 1 


n 


2 


2 1 


. 1 


K=.24+3.  25X-4.5X  Eqn  (42) 

x = h/b 

This  equation  enabled  a quantitative  analysis  of  induced 
drag  reduction  with  altitude  (for  a given  wing  span). 


D.  SUMMATION 

An  equation  for  power  was  optimized  in  a previous 
section  and  yielded  several  variable  terms,  of  these  HL/b 
was  the  dominant  term  and  an  equation  for  it  was  developed 
using  the  aircraft  geometry  and  Euler's  buckling  equation. 
A parasite  area  term  was  identified  and  an  expression 
developed  based  on  the  frontal  area.  Reductions  in  induced 
drag  were  noted  due  to  ground-effect  and  the  reduction  was 
quantified  according  to  a ratio  of  altitude  divided  by  wing 
span.  Oswald's  wing  efficiency  was  the  last  remaining 
variable  and  its  improvement  was  cited  as  possible  for  a 
fixed  sail  area  independent  of  the  WL/b  equation.  Design 
modifications  and  optimization  of  the  above  factors  were 
examined  to  bring  the  power  required  for  flight  within  the 
attainable  power  available  regime. 
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III.  SEARCH  FOB  OPTIMUMS 


A.  COMPUTATIONAL  ANALYSIS 

The  logical  approach  for  further  optimization  at  this 
point  appeared  to  be  siailar  to  the  optimization  technique 
of  the  power  equation  from  the  Optimum  Power  section.  Since 
the  intention  was  to  minimize  the  power-required,  it 
followed  that  further  optimization  of  the  power  equation  was 
the  obvious  selection  for  further  analysis  and  the 
development  of  an  expression  for  HL/b  was  developed 
accordingly.  The  abbreviated  terns  in  Equation  (36)  were 
expanded  and  it  was  intended  to  optimize  this  equation  with 
respect  to  span  Just  as  the  power  equation  had  been 
optimized  with  respect  to  velocity.  Setting  the  first 
denvat  : • • of  Equation  (38)  equal  to  zero  and  solving  for 
optimum  span  resulted  in  a tenth  order  equation  in  terms  of 
span.  Further  efforts  at  optimization  through  this 
technique  were  abandoned  and  an  extensive  analysis  of 
optimization  through  a digital  computer  approach  was 
selected  . 


B.  COMPUTER  ANALYSIS 

1.  Hi lot  £H  2£timizjtion 

The  first  computer  program  developed  was  called 
"Plot  P"  because  of  the  subroutine  employed  for  graphical 
output  of  HL/b.  It  was  felt  that  there  would  be  an  optimum 
value  of  HL/b  for  a given  set  of  inputs  as  the  span  was 
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allowed  to  vary.  The  required  inputs  to  this  program 

included  the  constants  K ,K  , h , T , and  sail  weight/unit 

d t t 

area.  Also  required  were  the  variable  terns  "Q"  (the 

inverse  of  the  inner  slenderness  ratio)  which  was  a lore 
convenient  numerical  expression  than  the  slenderness  ratio 
itself;  "X"  which  was  a ratio  of  the  outer/inner  spar  radii, 
and  "S"  the  sail  area.  Plots  of  BL/b  vs.  b produced  a 
niniaua  for  each  choice  of  the  paraneters  Q,X,and  S.  The 
analysis  then  becaae  a systematic  examination  of  these 
conbinations  for  an  overall  niniaua. 


2 . "Tab”  Outputs 

The  "Tab"  prograa  was  created  to  produce  a 

tabulation  of  the  paraneters  chosen  above  and  the  results 

that  were  obtained.  Input  data  included  the  sane 

inforaation  as  required  by  "Plot  P"  plus  the  estimated  body 

weight  (B  ) and  the  optiaua  wingspan  indicated  by  the  "Plot 
b 

P"  program.  The  critical  factors  in  the  "Tab"  program 

output  reduced  to  the  inner  and  outer  wing  spar  thicknesses, 
the  lift  coefficient  and  the  calculated  power  required. 
Other  output  provided  trend  comparisons,  spar  radii,  spar 
weights,  overall  weights  etc.  The  combination  of  these 
computer  programs  enabled  a systeaatic  examination  of  the 
variables  Q,X,  and  S for  the  best  possible  conbinations  to 
niniaize  Equation  (38)  while  remaining  within  the 
established  limitations. 

3.  Variation  of  Parameters 

Figure  12  illustrates  the  trend  in  BL/b  caused  by 
separately  increasing  each  of  the  parameters.  It  suggested 
improvements  in  BL/b  from  increasing  sail  area,  increasing 
spar  radii  and  from  the  outer  spar  radii 
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being  larger  than  the  inner  ones.  Increasing  spar  radii 

resulted  in  decreasing  spar  thicknesses  and  consequent 

weight  and  power  reductions.  Larger  outer  spars  were 

attributed  to  the  differences  in  respective  spar  end 

conditions  as  noted  in  the  wing-spar  coluan  analysis.  It 

seeaed  reasonable  that  a pin-ended  coluan  would  require  a 

larger  radius  than  a fixed-ended  coluan  to  successfully 

withstand  buckling  loads.  While  this  configuration  aight 

present  a slightly  unorthodox  appearance  it  was  considered 

desireable  froa  a design  standpoint  Pigure  13  shows  the 

paraaeter  influences  in  broader  teras.  Data  for  both 

Figures  12  and  13  were  obtained  froa  the  "Tab"  and  "Plot  P" 

coaputer  prograas.  Once  the  wing  spar  thickness  had  been 

chosen,  the  Q and  X paraaeters  were  adjusted  to  attain 

coaputed  spar  thicknesses  as  near  to  .030"  as  possible  in 

both  inner  and  outer  spars.  This  procedure  was  repeated  for 

sail  areas  ranging  from  500  to  1500  square  feet.  Spar 

thicknesses  equal  to  ,030"±.001"  under  various  conditions. 

Wing  spars  were  the  optiaua  spans  indicated  froa  "Plot  P" 

coaputer  output.  WL/b,t  ,t  ,Cl,  and  power  (including  ground 

i o 

effects)  were  calculated  by  "Tab."  Having  satisfied  the 
spar  thickness  criteria,  the  lift  coefficients  and  power 
required  values  were  examined.  All  power  values  appeared  to 

WL/b  P b t t 

g i o 

<5  » l t l l 

to  1 1 l t t 

U 4 J t t I 

Figure  13  EFFECTS  OF  SEPARATE  VARIATION  OF  PARAMETERS 
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be  within  a feasable  range  as  did  most  of  the  lift 
coefficients.  It  appeared  at  this  point  that  further 
improvement  in  HL/b  could  only  be  achieved  through  an 
increase  in  Oswald's  efficiency  factor,  a decrease  in  the 
parasite  drag  or  a reduction  in  the  body  weight.  The  latter 
had  been  only  a first  approximation  and  not  yet  the  subject 
of  detailed  examination.  Fuselage  spar  weights  were 
calculated  using  a third  computer  program  which  also 
calculated  weight  and  balance  information,  as  well  as  shear 
and  bending  stresses  in  the  fuselage  spar.  Because  the 
program  also  dealt  with  different  sail  shapes  the  entire 
discussion  was  treated  in  a subsequent  section. 

C.  TRADE-OFF  ANALYSIS 

Several  options  remained  available  at  this  point  which 
might  serve  to  increase  the  power  available  (2-man  crew) , 
decrease  the  aircraft  weight  (lighter  construction 
materials)  or  increase  the  wing  efficiency  (wing  shapes 
other  than  pointed  wingtips) . In  addition,  the  parasite 
drag  and  sail  weight  sensitivities  were  examined.  Computer 
output  lists  comparisons  of  various  trade-off 
considerations. 

1 . Two-Ban  £r ew 

The  idea  of  a two-man  crew  was  entertained  briefly. 
Its  attractive  feature  was  a doubling  of  the  power  available 
to  propel  the  craft.  As  indicated  by  the  optimum  power 
equation,  however,  the  relationship  between  power-required 
and  weight  was  not  a linear  one.  The  geometric  penalties 
incurred  from  another  crewman  more  than  offset  the  linear 
increase  in  power  available.  Thus  a single-piloted  craft 
was  selected  relatively  early  in  the  study  and  the 
multi-piloted  version  was  not  considered  further. 
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2.  Magnesium  versus  Aluminium 


Magnesium  was  considered  as  a structural  alternative 
to  aluminium  principally  because  of  its  significant  decrease 
in  density.  At  the  sane  time  there  was  also  a decrease  in 

the  Modulus  of  Elasticity  (a  measure  of  a material's 
stiffness)  which  might  adversely  affect  the  structural 

I 

integrity  of  the  craft.  Also  allowable  shear  and  bending 
stresses  were  lower.  "Plot  P","Tab"  and  "Mt-Balance" 
programs  were  run  using  magnesium's  properties  and  the 
results  compared  to  aluminium.  Equivalent  spar  thicknesses 
in  both  metals  required  larger  magnesium  spars  (lower  Q and 
larger  X)  and  slightly  longer  optimum  wing  spans  were 
indicated.  It  appeared  that  a decreased  modulus  of 

elasticity  dictated  an  overall  increase  in  spar  size. 
Improvements  in  magnesium  spar  weights  were  offset  by  the 
larger  and  heavier  fuselage  spar.  The  shear  and  bending 
stresses  appeared  to  be  critical  factors  in  magnesium.  An 
interesting  combination  might  include  an  aluminium  fuselage 
with  magnesium  wing  spars.  Further  investigation  of  the 
magnesium  option  is  warranted  concerning  maximum  shear  and 
bending  stresses. 

Aluminium  was  considered  the  superior  selection 
based  on  the  above  analysis.  It  had  the  additional 
advantages  of  availability  and  lower  costs. 

3.  gacaaiia.  2 tag  sansiiixiii 

In  order  to  assess  parasite  drag  influences  on  the 

power  required,  the  value  of  "f  " was  set  to  zero.  This 

man 

drag  modification  was  introduced  into  the  "Tab"  computer 
program  and  its  results  compared  with  a run  in  which 
2 

f =3.5  ft  as  originally  assumed.  For  the  1000  foot  sail 
man 

area  the  power  required  dropped  from  .345  to  .325 
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horsepower.  In  otherwords  a complete  elimination  of  the 


f drag  tern  yielded  a reduction  in  power-required  of 
■an 

approximately  6%.  This  confirmed  earlier  suspicions 
concerning  the  drag  variable  as  a fairly  insensitive  factor. 
Efforts  directed  at  enclosing  the  pilot  in  a streamlined 

fuselage  would  not  eliminate  f and  could  be  costly  in 

■an 

terms  of  increased  weight. 

4.  Sail  Heigh*  Sensitivity 

Original  conceptions  for  sail  weights  consisted  of 
the  lightest  dacron  material  available.  Dacron  was 
considered  preferrable  to  nylon  because  of  the  former's 
resistance  to  stretch  under  loaded  conditions.  A weight 
figure  of  3/4  ounce  per  square  yard  was  provided  by 
commercial  sources  and  was  utilized  in  early  sail  weight 
calculations.  Later,  information  from  Sailrite  Kits,  a firm 
dealing  in  dc-it-y ourself  sails,  indicated  that  2.2  ounce 
dacron  might  be  the  lightest  obtainable  and  that  this  figure 
was  for  a 36"x  28.5"  section.  This  revised  sail  weight  was 
introduced  into  the  "Tab"  computer  program  to  examine  its 
effect  on  power  requirements.  Not  surprisingly  the  greatest 
deviations  were  noted  in  the  larger  sail  areas  with  the 
range  of  deviations  running  from  a low  of  about  a 4%  power 
required  increase  to  10%  on  the  1500  square  foot  sail.  The 
1000  square  foot  sail  required  about  an  8%  increase  (.317  vs 
.345  hp) . Thus,  a rather  dramatic  increase  in  the  weight  of 
the  sail  material  was  not  catastrophic  to  the  overall 
performance  requirements.  It  did  however,  indicate  that 
lighter  sails  would  be  strongly  preferrable  to  the  2.2  ounce 
case.  Further,  while  the  1000  square  foot  sail  area  was 
clearly  superior  for  the  lighter  sail,  a 750  square  foot 
sail  would  be  equally  effective  in  the  heavier  instance. 


The  latter  assumed  that  the  higher  lift  coefficient  was 
acceptable. 

5.  Taper-Batio  versus  Pointed  Wingtip 


Oswald's  efficiency  includes  the  products  of  the 
planform  efficiency  (which  is  1.0  for  elliptical  wings), 
body  interference  effects,  separation  and  stall  effects, 
tail  contributions,  propeller  effects,  flap  effects  etc. 
Experience  has  shown  that  an  Oswald  efficiency  of  85%  is 
about  the  best  that  an  elliptical  wing  can  hope  to  obtain. 
Further  decreases  occur  as  wing  shape  is  modified.  In 
practice  elliptical  wing  construction  is  a difficult  matter 
and  truly  elliptical  wings  are  rarely  built.  Approximations 
to  an  elliptical  wing  can  be  obtained  by  a wing  with  a 
finite  wingtip  chord.  Glauert  Ref  ' 3”  has  shown  the  effects 
of  wing  taper  ratio  (wingtip  chord  divided  by  wing  root 
chord)  on  fractional  drag  increase  over  an  elliptical  wing 


(Figure  14)  From  aerodynamics  Cd*C^  /irARe  and  therefore 
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Eqn  (43) 


Figure  14.  TAPER  RATIO  EFFECTS  ON  DRAG  INCREASE 

It  was  apparent  that  some  benefit  could  be  gained  from  using 
a tapered  airfoil  over  a pointed  wingtip.  The  question 


posed  a tradeoff  between  a weight  increase  froa  the 
additional  wingtip  booa  requirements  versus  anticipated 
gains  in  wing  efficiency.  The  wingtip  booa  would  have  to 
sustain  bending  aoaents  froa  the  lifting  force  of  the  sails. 
Hith  this  consideration  in  aind  a taper  ratio  of  20%  was 
selected  as  still  having  a favorable  reduction  in  drag  while 
also  providing  a saall  tip  chord,  oaaller  tip  chords  would 
decrease  the  bending  aoaents  with  consequent  weight  savings 
in  the  wingtip  boon.  A coaputer  prograa  was  created  to 
exaaine  this  question  and  its  daacri  ■♦ion  follows. 

a.  "Height-Balance”  Coaput*-  Prograa 

This  prograa  was  initially  created  to  solve  the 
weight-and-balance  question  relating  ' ■ filot  and  propeller 
positioning.  It  was  expanded  to  calculate  a fuselage  weight 
based  on  shear  and  bending  stresses  and  further  expanded  to 
exaaine  a wingtip  booa  in  a siailar  Banner.  Inputs  required 
froa  the  "Tab"  coaputer  prograa  included  the  coaputed  spar 
weights  and  the  wing-lift,  and  froa  "Plot  P"  the  optiaun 
span.  In  addition  the  saae  assumptions  were  required  plus 
any  desired  wing  taper  ratio.  Outputs  included  the 
positioning  of  the  aan  and  propeller,  size  and  weight  of  the 
fuselage  spar,  sail  diaensions,  an  accurate  body  weight  and 
the  size  and  weight  of  the  wingtip  booa.  Shear  and  bending 
stresses  were  also  listed  as  a aatter  of  interest.  Naxiaua 
shear  and  bending  values  were  obtained  froa  Ref  '5"  using  a 
yield  stress  for  a hollow  tube  (good  for  thicknesses  up  to 
.249  inches).  Ref  '5"  further  indicated  that  "compressive 
yield  stress  was  approxiaately  the  saae  as  tensile  yield 
stress"  for  this  case. 

Optimum  parameter  cases  with  sail  areas  ranging 
froa  500-1500  square  feet  were  examined  with  detailed 
analysis  on  the  1000  square  foot  sail.  In  this  instance 
arbitrarily  selected  fuselage  radii  of  2"  and  3"  were 
examined  with  each  of  these  considered  for  tapered  and 
pointed  wingtips.  The  smaller  fuselage  spar  radius  was 
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found  to  is  adequate  in  terms  of  the  stress  analysis  and 
preferrable  in  terms  of  weight  savings.  Tapered  wings 
produced  lighter  and  smaller  fuselage  spars  because  the  root 
chord  was  reduced  and  many  assumptions  on  fuselage  length 
were  in  terms  of  the  root  chord.  A weight  penalty  of  three 
pounds  for  both  wingtip  booms  for  the  1000  square  foot  sail 
appeared  to  be  well  within  acceptable  weight  penalty  limits. 
A rough  estimate  of  drag  decrease  was  obtained  from  Figure  9 
comparing  a pointed  wingtip  with  a 20%  taper  ratio  and  both 
of  these  to  an  elliptical  wing.  The  curve  peaked  on  the  left 
side  of  the  curve  at  about  13%  . This  was  subtracted  from 
the  85%  efficiency  of  the  elliptical  wing  and  for  pointed 
wingtips  an  efficiency  of  70%  was  used  in  calculations.  For 
the  tapered  wingtip  estimates  yielded 
(1 . 13-1 . 04) /I . 13=. 0796,  or  roughly  an  8%  drag  reduction  over 
the  pointed  wingtip.  Then  .0796  times  .08  yielded  about  5% 
improvement.  Values  of  .75  were  used  for  the  tapered  wings. 
Since  induced  drag  and  wing  efficiency  were  seen  to  be 
inversely  proportional  from  Equation  (43) , then  a 5% 
reduction  in  induced  drag  might  suggest  a 5%  increase  in 
wing  efficiency.  Improved  wing  efficiencies  were  introduced 
into  "Tab"  to  examine  the  effects  of  improved  wing 
efficiency  versus  a wingtip  boom  weight  addition.  Results 
indicated  an  decrease  in  power-required  of  approximately  6% 
for  the  1000  square  foot  sail  and  slightly  less  improvement 
in  others.  It  was  therefor  concluded  that  the  tapered  wing 
was  a desireable  option  due  to  the  decreased  power 
requirements  despite  the  increase  in  weight  and  structural 
complexity. 


0.  AREAS  FOR  FORT  HER  CONSIDERATION 

While  the  structural  design  of  the  craft  had  been 
determined  at  this  point  several  areas  worthy  of  further 
consideration  exist  and  exploration  would  be  desireable,  if 


not  essential,  prior  to  actual  flight  testing  such  a craft. 

1 • Propeller  Optimization 

Given  the  huaan  capacity  for  sustained  power  output 
an  analysis  of  the  optimum  rpm  and  propeller  blade  diameter 
would  enable  an  increase  of  the  power  efficiency.  The  power 
efficiency  for  all  calculations  thus  far  was  chosen  to  be  .8 
. Naturally  any  improvement  in  this  would  further  decrease 
power  reguirements.  Reference  '8"  discusses  this  problem  in 
depth.  Due  to  the  design  of  the  aircraft  to  this  point  a 
propeller  diameter  of  approximately  eight  feet  appeared  to 
be  about  a maximum.  Prom  Sherwin's  *8"  figure  reproduced 
below  an  rPM  of  about  250  would  yield  a propeller  efficiency 
of  better  than  8554.  This  would  be  a significant  addition  to 
the  improvement  of  the  propulsion  efficiency. 


.<10 


.30 


RfVi 
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V i 17 

Figure  15 


PROPELLER  EFFECTIVENESS  WITH  VARIATION  OF  LENGTH  AND  RPM 
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2.  Stability  and  Control 

V 

Little  attention  had  been  accorded  the  actual  control 
problem  of  the  aircraft  and  its  gust  response  would  have  to 
be  analysed  in  some  detail.  Plutter  problems  would  need  to 
be  addressed  as  well,  since  the  aircraft  was  essentially 
two  large  moments  of  inertia  (wing  and  horizontal  tail) 
there  remained  a question  of  torsion  due  to  gusts  and 
asymmetrical  lifting  forces  from  the  wings.  If  a wake 
disturbance  from  the  wing  arrives  at  the  tail  in  a time 
comparable  to  the  time  for  one  cycle  of  a torsional 
oscillation  between  wing  and  tail  then  flutter  instability 
can  develops.  Purther  restraining  stays  would  be  essential 
to  eliminate  this  problem.  Finally  a small  prototype  should 
be  constructed  to  examine  the  validity  of  all  assumptions 
and  the  control  responses  necessary.  In  addition, 
uncertainties  as  to  the  sail's  actual  performance  remain  to 
be  resolved. 


IV.  CONCLUSIONS 


The  concept  of  a sail-wing  and  exposed  structure 
appeared  to  produce  a vehicle  possessing  performance 
characteristics  competitive  with  conventional  rigid-wing 
designs.  However  , while  the  optimum  power  equation 
indicated  areas  for  optimization  analysis,  man-powered 
flight  is  impossible  without  ground  effects. 

For  each  sail  area  selected  an  optimum  combination  of 
wing  and  fuselage  size  and  shape  existed.  Principal 
improvements  were  seen  to  be  in  the  area  of  weight  reduction 
although  other  refinements  were  beneficial.  It  appeared 
that  a light  sail,  tapered  wing  and  magnesium  wing  spars 
with  an  aluminium  fuselage  spar  would  give  the  highest 
performance  for  a sail-wing  type  of  aircraft. 

While  many  details  of  construction  engineering  have  not 
been  examined,  indications  are  that  the  aircraft  flying  at 
Kitty  Hawk  might  ave  been  closer  to  a successful  man-powered 
prototype  than  anything  thit  has  yet  followed. 
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